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By using a Path Integral formalism, we study the fluctuating fields in the order parameter 
for a 3D uniform Fermi gas at T < Tc. The BCS-BEC crossover condition is included into 
the analysis in the framework of the usual s-wave scattering approximation. The general 
expressions for the temperature dependence of the partition function and the quadratic 
average in the fluctuations fields are explicitly established in terms of the BCS-Gor'kov 
correlation propagators. 
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2. Introduction 

The problem of the crossover from BCS theory with cooperative Cooper pairing to the for- 
mation and condensation of composite bosons has attracted considerable attention after 
the discovery of the high T c superconductors. Monte Carlo simulations for the attractive 
(negative U) Hubbard model have shown that, in the intermediate coupling, the crossover 
region displays highly anomalous correlations in degenerate Fermi system above T c pQ. 
More generally, the crossover phenomenon is also of interest for other problems: excitonic 
condensates, superconductor-insulator transition and Cooper pair formation in resonant 
Feshbach channels |2],|S]- 

By using the path integral formalism, we study a 3D continuum model of fermions with 
attractive interactions, at finite temperature. In the first part, we perform a step-by-step 
derivation of the Time-Dependent Ginzburg Landau (TDGL) equation for the particular 
case in which the order parameter vanishes (T = T c ), and arbitrary strength coupling g. 
We show that in the weak coupling limit, the coefficients in TDGL equation will reduce to 
the classical ones obtained by Gor'kov (1959) in the framework of the microscopic approx- 
imation. This formalism is extended for the case T < T c , the s-wave pairing symmetry 
and the Mean-Field Approximation (MFA), following closely the approach of Randeira's 
group, reported in Ref. 0j. 

The main result in this report is the explicit calculation of the quadratic average fluctu- 
ation field (| 7] | 2 ), where 77 stands for the fluctuation from the mean field value in the 
order parameter A . We obtain the complete set of self-consistent equations for the dy- 
namical (oj 7^ 0) case for T < T c , and one expression for (| rj | 2 ), by keeping only the zero 
frequency component (the static approximation). Further developments are proposed in 
terms of numerical analysis, high-order expansion (terms up to fourth order in the effec- 
tive action), modification of the "starting" Hamiltonian (i.e. by adding the terms that 
describe both the Feshbach resonance effects and preformed Cooper pairs). 
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3. The Hubbard-Stratonovich Transformation 

The Hubbard-Stratonovich (HS) transformation maps certain classes of interacting fermion 
problem, onto non-interacting electrons moving through an effective field. Let us consider 
that the interaction part of the Hamiltonian density involves the product of two bilinear 
terms: 

Hi = -gV, (x) ^ T (x) ^ (x) ^ (x) ^ T (x) . (1) 

Taking into consideration the bilinear form for the pair density A (x) = (x) (x), the 
Hamiltonian for the interacting part in the continuum space is given by: 



H T 



or 



J Hd 3 x = J dxA^ (x) A (x) , (2) 
ff/ = -</£4 4 J ( 3 ) 



3 

on a lattice. 

The HS transformation replaces the combination — gA (x) — > A (x), in the following way: 

-gA^A -> A^A (x) + A(x)A (x) + A ^ . (4) 

The field A can be though as the exchange boson which mediates the original interaction. 
In this point is necessary to recall several useful identities of complex integration. Bose 
coherent states are parametrized by complex vectors z = (-21,-22, • • ',zn), z % — x i + iyi, 
where are real and N is the number of single-particle states. The complex integration 
is defined as: 

d 2 z = j 00 n (^^) • (5) 

A useful identity of complex integration over a single variable z is: 

1 1 rf 2 ^ (z*) n ^ m exp = < m . (6) 

Using this definition, it follows that for any complex matrix G, whose Hermitian part has 
only positive eigenvalues, the Gaussian integral is given by: 

J d 2 zexp (-z*Gz - z* a z - z*z b ) = det | G | _1 exp (z* a G~ 1 z b ), (7) 
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where z* a , z& are any complex vectors. The formalism for fermionic fields is carry out in 
a similar way (this time the Zj compounds are Grassman numbers), and Eq. (7) can be 
unified by the single expression: 

J d 2 z exp (-z*Gz - z* a z - z*z b ) = det | G | _J? exp (z k a G~ 1 z b ) , (8) 

where rj = 1 for bosons and rj = — 1 for fermions. 

Under the definitions A = Ai — zA 2 , A = Ai + zA 2 , dAdA = 2idAidA 2 and with the 
identities: 

/ dA^A 2 e(- A ?- A ^ = ng, [ ^V AA /* = 1, (9) 

J J Zntg 

where A (x, r) = A (x) is a function of space and time, defined on a grid in space-time 
that can be made infinitely fine. The natural generalization of Eq. (9) can be written as: 

J T)[A,A}e-f d3 *f° dTA ~ (x)A{x)/9 = 1, (10) 

where f3 = 1/T and r represents the imaginary time in the finite temperature formalism 
The integration measurement X>[A, A] is defined as: 

D[A,A]^n rfA(Xj,r ^ A(Xj,r) > (H) 

where Af = 2nig/ (A 3 xAt) is a normalization and we have taken the continuum limit 
(Ax, At) — > 0, replacing the discrete sum over the grid r . A (xj)A (xj) by the inte- 
gral / c/ 3 x J'q drA (x)A (x). 

3.1 Constructing the generating integral 



The imaginary time generating functional is defined as: 



Z = TrT r exp 



- / dr I d 3 x.H 



o 



(12) 



where Tr{-} represents the Trace operator over the quantity (■), and TC is the Hamiltonian 
density of the system. By manipulating the Tr and T T operators, the partition function 
can be written as jjj: 

Z = Jv^^je- 8 ^^), (13) 
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with: 



S ({*}{*})= J I d 3 x[y a (x)d T V a (x)+H], 



as the action functional in the fields and 



^ CT (x) + Hi, 



(14) 



(15) 



where Tii corresponds to the interacting part for the Hamiltonian density defined in Eq. 
(1), and (j, is the chemical potential that fixes the average density n. 
Introducing the unit identity in Z: 

Z = J Z>[tf,tt] j V[A, Aje-^^H*})-^^^*^*). (i 6 ) 

By inserting the explicit form of S , the exponential argument in Eq. (16) reads: 

S (WW) + - / rf 3 x I drAA = 
9 J Jo 



dr 



V 2 \ - - 1 - 



Under the transformations: 

A (x) -> A (x) + gA (x) ; A (x) -> A (x) + ^ (x) , 
the last two terms in Eq. (17) become into: 

-gA ] A + - AA -> -AA + AA + A f A = -AA + A^, ^ T + . A. 
9 9 9 

This construction must be consistent with the motion equations: 



2m 



gVf&pV,, = 0, 



(17) 



(18) 



(19) 



-d T ~ — ~ 11 ) ~ g^a^^p = 0, 



(20) 
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derived from the conditions SS/S^ a = and 5S , /5\I / cr = |T2|. Then, the l.h.s. of Eq. 
(17) can be finally written as 1 : 



eTx / dr 



Iaa 

9 



(21) 



In this procedure we have absorbed the interacting term by replacing it into an effective 
field which couples the bilinear structure A. Integrating out over the fermionic fields, the 
partition function takes the form: (Cfr. Eqs. (8) and (16)) 



Z= Z>[A, A]e" 



M 



:22i 



where S e fj is the effective action which can be written in terms of the (inverse) Nambu 
propagator G -1 (defined below). The intermediate steps are developed in the following 
way through the redefinition of Eq. (21): 



with 



Recalling Eq. (16): 



Z 



G 1 



A 



A 



(23) 



(24) 



P[A,A] 



(25) 



By using Eq. (8) and the expression det (M)=e^ rlnM , the effective action can be rewritten 

as: ^ 

S eff [A,A] = [ d 3 x [ dr [-AA - Tr In G 1 (A)] . (26) 
J Jo 9 

4. Gap Equation Derivation and cut-off criterion 



The last integral is analized in the context of static, uniform-saddle points and fluc- 
tuations about them. The trivial saddle-point A = 0, which is stable at sufficiently high 
T for all couplings, will become unstable below a certain temperature denoted as T . The 

1 In this step is important to notice: ^ , cr 3 r v I' a = ^ a d T ^ a . 
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transition temperature is defined by the solution of SS e ff/SA |a=o— 0, or, in equivalent 
way: 0. M- 

- = ^tanh(£ k /2T )/2£ k . (27) 



9 k 



:28i 



Derivation of Eq. (27) 

By applying spatial Fourier transformation in the second term of Eq. (26), we have: 

S eff [A, A] = J* dr (X A A - In HJ det [-d T - h kjT ] J J , 

and h k T is a matrix which contains arrays of kinetic energy terms £ k = £k — with 
= k 2 /2m, as well as the coupling field A: 

Kr=(*k if)- (»> 



-A 

The time- Fourier transformation \l/ k (r) = /3~ 1//2 ^k in ^~ lulnT \ provides: 

S eff [A, A] w^AA-lnJIdet^-hJ, (30) 

" k,n 

where the imaginary-time fluctuations dependence in A have been neglected, and {uj n } 
represents the set of fermionic Matsubara frequencies defined by: u n = (2n + 1) 7rT, with 
n as an integer number. Further manipulations gives: {hi]Xi(' ' ') = Xm m (' ' ') anc ^ 
| det [iu n - h k ] \=ul + £l + A 2 } 

Seff [A, A] = ^AA - ^ K + I A | 2 ). (31) 

" k,n 



Using the condition 5S e ff/5A = 0: 

5S eff [A,A] _ VP A 



k,n 

or: 



i i ^ i 



E ,„ 2 ^^IA| 2 =^ (32) 



1 1 V 1 

The complex integration over closed contours give us a simple formula for evaluating sums 
like (33) |JT] 

— F {iuo n ) = sum of residues of F (oS) f (uj) at poles of F (u), (34) 
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where / (uj) = (l + e^ w ) 1 is the Fermi distribution. Clearly, the poles in F (u) are ±E^, 
with Ek = | A | 2 , therefore the expression (33) can be rewritten as: 

vL lR 9*T~ ■ (35) 



g V ^ V 2^ k 2£ k 

Straightforward calculation leads to: 

1 1 ^tanh(/3£ k /2) 

rv? 2E k • (36) 

which is reduced to Eq. (27) when A = 0. 

A BSC type cutoff cannot be used to access the strong coupling Bose regime. The use 
of the two-body s— wave scattering theory regulates the ultraviolet divergence in the gap 
equation through the relation: 



rn 



4vra s g ^ 2e_ 

where a s is the scattering length, and e k = k 2 /2m. In the following section, we will derive 
the regularization condition (37) from the elementary s-wave scattering theory. 

4-1 Derivation of the Regularization condition in the Effective Interaction 

The concept of the scattering amplitude involves only the directions of the initial and 
final momenta of the particle undergoing scattering ^2] • It can therefore be naturally ap- 
proached by formulating the scattering problem in the momentum representation, where 
there is no question of the spatial distribution of the process. 

Let us transform to the momentum representation the original Schrodinger equation 

-|^V 2 ^ (r) + [U (r) — E]ip (r) = 0, (38) 
changing from coordinate to momentum wave functions, i.e. to the Fourier components: 

a (q) = J if, (r) e~ iq W. (39) 

Conversely, 

1> W = / -^U (q) e iq r . (40) 
J (2vr) 
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By multiplying (38) by e iq r and integrate over dV. In the first term: 

d 3 q' 



J dVe~ ici r X7 2 4> (r) = J dV J e~ iqr x J 



a(q')(-9 / ) 2 e iq '- r 



/^H ,2 )«(q / )5 (3) (q'-q) = -A(q) 



(41) 



The substitution of Eq. (40) into the second term of Eq. (38) 1 gives: 

d 3 q' 



U (r) V (r) e" jq T dV 



(2k) 



^(q-q')a(q'). 



(42) 



Thus Schrodinger's equation in the momentum representation becomes: 

!£- £ >W|J f/(q - q ' )a(q,)=o ' 



(43) 



The wave function describing the scattering of particles with momentum Tik has the form: 

^ (r) = e* kr + Xk (r) , (44) 

where Xk (r) is a function whose asymptotic form (as r — > oo) is that of an outgoing 
spherical wave. Its Fourier component is: (see Eq. 39) 



«k (q) = 5 (3) (q - k) + X k (q) 

and substitution in (43) gives the following equation for the function Xk (q): 
{{q 2 - A; 2 )5 3 (q-k) = 0) 



(45) 



2m 



{k 2 -q 2 ) Xk(q) = ^(q-k) + J (q - q') X k (q') , 



(46) 



where E = h 2 k 2 /2m. This equation may conveniently be transformed by using instead of 
Xk (q) another unknown function defined by: 



Xk(q) 



2m F (k, q) 
Wq 2 -k 2 -i0' 



(47) 



iHere, we must insert U (r) = / -0^U (k) e lk r , and (2tt) 3 <5< 3 ) (k + q' - q) = / d^e l ( k+q '" q ) . 
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This eliminates the singularity at q 2 = k 2 in the coefficients of (46), which becomes: 

d 3 q' U (q - q') F (k, q') 



F(k,q) = -[/(q-k)-|? J 



(2tt) 3 q' 2 -k 2 -i0 



(48) 



The term iO, which denotes the limit of i8 as 5 — > +0, is included in the (47) definition 
in order to give a definite sense to the integral (48), since it establishes the manner of 
passage round the pole q' 2 = k 2 . The Fourier transformation gives: 



2m r d 3 q F (k, q) e iq r 
~W J (2tt) 3 q 2 - k 2 - 20 



xkW=^ : 8 ; (49) 



Using <i 3 g = q 2 dqdQ n , and first integrate over <if2 q i.e., over the directions of the vector q 
relative to r: 

Xk (r) = " 2rfV /_ g»-*»-i0 - (50) 

where n' = r/r. The integrand has poles at points q = k + iO and q = —k — iO; the path 
of integration in the complex q— plane passes respectively below and above these values. 
The integral tends to zero as r — > oo, and the integral round the loop is given by 2ni 
times the residue of the integrand at the pole q = k. The final result is: 

xM = ^—F{knM), (51) 

where n is a unit vector in the direction of k. We have derived the requiered asymptotic 
form of the wave function, and the scattering amplitude is: 

Tfl 

f(n,n') = ^F(kn,kn'). (52) 

In the first approximation, omitting the integral term in (48), we have F (k, q) = — U (q — k). 
In the next approximation, we substitute it in the integral equation, then: 



/(n,n') = — ^ 



2m f d 3 k" U (k' - k") U (k" - k) 



(2tt) 3 k 2 - k" 2 + iO 



(53) 



with k = kn, k' = kn'. The first term is often called the Born approximation. In the 
second approximation, the scattering amplitude does not have the symmetry property 
/ (k, k') = /*(k', k). At first sight it may seem that the integral term in (53) is also 
symmetrical with respect to interchange of the initial and final states. However, this 
symmetry does not exist, because the path of integration and the direction of the passage 
round the pole are altered when the complex conjugate expression is taken. 
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In order to make these ideas quantitative, consider the the two-particles scattering, where 
the particles are assume to have equal masses m, and therefore m T — m/2. In this case, 
the scattering matrix F (k, q) in Eq. (48) satisfies the so called Lippmann-Schwinger 
equation [13 : 

F (k, q) = -U (q - k) - / ^^-^^^ (54) 
V ,HJ VM J J (2vr) 3 hq' 2 /m- E - iO ' V ; 

where E = h 2 k 2 /m is the energy eigenvalue and U (q — k) is the Fourier transform of the 
bare atom-atom interaction. The scattered wave at large distances and for zero energy 
(E — 0, k — 0), may be calculated from (54). In this case, the previous formulation must 
be slightly modified through the change 2m — * m [14J. Using the formula: 



fl3U pik(r-r') i 

^> (55) 



Xo(r)«-^F = -^. (56) 



(2tt) 3 k 2 4tt 
expression (49) can be approximate by: 

mU(0) __ a 
A-nh 2 r r 

In the context of Born approximation, which is obtained by taking only the first term on 
the right side of (54), the scattering length is given by 

= ~4rf~' (57) 

The relation (57) between U and the scattering amplitude is not exact, but is only valid 
up to first-order terms. We are also interested in higher order terms in the energy, the 
formula (57) has to be corrected. Using Eq. (54): (F w —Airha s /m) in the l.h.s, and 
F w -U (0) in r.h.s) 

Anh 2 a s U 2 (0)^ 1 



m V ' hq 2 /m 

q 

Further arrangement finally leads to the expression (37): 

m ~ 1 | 1 y 1 m 

47ih 2 \a s \ V^Y V ^nq 2 jm 1 ] 

Eliminating the coupling constant g = U (0) in Eqs. (27) and (59), we obtain: 

m v / tanh(^ q /2To) 1 



(60) 
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In order to solve the last equation for T , we must determine the chemical potential [i as 
a function of the coupling interaction g and the temperature, by using the relationship 
N = —dQ/dfi. The saddle point approximation for the thermodynamic potential, Qq = 
S e jf [A] /P, leads to the number equation: (see Eq. (31)) 



n <— (T ' = ~dH = -p g c^+IAI 2 = ? ^ tanh IT J ' (61) 
with E k = a/^ + A 2 - F° r bosons: 

„(T^)=X;[l-^tanh^)]. (62) 
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5. Expansions in the Action Functional and TDGL Equation 



In order to obtain the contributions on the fourth-order fluctuations in the effective ac- 
tion S e ff [A, A] due to the bosonic fields A (A), we must start by writting the Nambu 
propagator (24) as: [T3] 



G 



-d T - JCa z + Act + + Act - = G 



-i 



A 



a 



A 



a 



(63) 



where K, = — V 2 /2m — /x, cr ± = (a x ±icr y ) /2, where <x,'s are the Pauli's matrices, and 
Go -1 is the non-interacting propagator. The expression for S e ff in Eq. (26) becomes: 



S eff [A, A] = J d 3 xjf dr 



A 



- Tr lnGfT 1 - Trln [l + G (Act + + Act - )] 



(64) 



Using the identity Tr In (1 + A) — — X^°=i (— 1)" Tr (A) n /n, we obtain the combination: 
Trln [1 + Go (Act+ + Aa~)] -> — Tr (G a+G' a- AA' + Goct~G ct + A'A) , (65) 

for the second order expansion in the last term of (64) (o" ± Gocr ± = 0). Rewritting (64) in 
terms of the Fourier variables {ioJ m , iu n , k, q), with iuj m = 2nT (m + 1/2), iiv n = 2irTn) 
as the Fermionic (Bosonic) Matsubara frequencies, the Green's function reads: 



Go 





{iu m + £ k 



(66) 



and the effective action becomes into: 
S eff [A, A] =S +J2 



A (a) I 2 1 

-A^+-Tr (G 0fc a + G _ fe+ ^- + G _ fc+g a~ G o^ + ) \ A (q) | 2 
9 ^ 

(67) 

J d 3 x J Q (irTrlnGo -1 , and the Tr operator involves the internal-vertex 
sumation over the fermionic four-momenta k. The operation | Tr (• • •) finally gives: 

' ^ " - - (6g) 



with Sn 



-Tr i 



^ (iu m - £k) 1 (iu m - iv n + £_ k+q ) 1 



Expression (67) can be reduced to: 

S eff [A, A] =S + J2 



r(g) 



(69) 
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for the Gaussian approximation. The coefficient T is defined by: 



-1-1 



m 



Ana. 



+ ^ \( iuJ m ~ fk) 1 (iUm - IVn + £-k+q) 1 + — 



(70) 



where regularization condition for the effective s-wave interaction g (See Eq. 37) has 
been inserted in the last expression. The resulting equation for the thermodynamic po- 
tential Q = Q — jS^ 1 J2 q l n T^ 1 is identical to the diagrammatic result of Nozieres and 
Schmitt-Rink (NSR), in the framework of particle-particle scattering t-matrix theory |16j . 

5.1 Static Part 

Expansions for the static compound of (q, 0) can be written as: T _1 (q, 0) ~ a + 
c | q | 2 /2m + • • •, where a is evidently given by: 



a = T' 1 (0,0) 



m 



Ana* 



^ (iw m - Ck) 1 (iu m + £-k) 1 



2e, 



m 



V [— - — 



X = tanh 



06 



(71) 



The calculation for the c-coefficient starts from Eq. (70): 



:q,o) 



m 



47ra> 



[ J- 1 -n F (&) - n F (C-k+q) 



L2e k 



Ck + £-k+q 



(72) 



Using the relations: 



1 - tip (£ k ) - n F (C-k+q) 



tanh ( j + tanh ^ (£ k + /i) 



/i = -^(g 2 -2g(k.Ji)), 
2m 

valid in the small momentum q approach, Eq. (72) becomes into: 



l-n F (£ k ) - n F (f-k+q) 



2 A" + -/3/iY - C^Ph\ XY 



(73) 
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where Y = sech 2 (/3£k/2). It is also useful to consider: 



h h 2 
+ 



Straightforward algebra leads to the O (q 2 ) coefficient in T _1 (q, 0): 



2m 



X (k-n) 2 X pY [3(k-n) 2 Y (3 2 (k-n) 2 XY 



2m 



5.2 Non-Linear term 



(74) 



(75) 



The fourth order fluctuations in O (A 4 ) for the effective action S e ff can be calculated 
in analogous way. In this case, we must analyze the non-zero contributions: 



G <7 + Go<7 Go<7 + Go"<7 + permutations 



<?2 <?3 
Gk+Q 2 



G-k+Qi 




G. 



k+Q 3 



Qi 



G 



k+Q 4 



94 



In the our case, we have interest in the q\ = q2 = q% = coefficient: b (0, 0, 0) which is 
defined by: 

— -Tr (GofeO- + Go-fc (T Gojt^Go^o" + Go_ k o~ Qo k o- + Qo_ k o Go k o~ + ) = 



-\ ( iuj m - &)~ 2 (ju m + Ck)" 2 = -\b (0, 0, 0) . 



(76) 
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The expression (69) must be formally rewritten as: 



S eff [A, A] =^o + E ' r(g) |2 + ^ £ MO, 0, 0) A (ft) A ( fe ) A (ft) A ( ft + ft - fe) . 

(77) 

X 



91,92,93 



with 



6 = 6(0, 0,0) = J] 



4er 



7f 



To study the evolution of the TDGL we recall the relationship 5S e ff/SA* = near to T c 
and for A slowly varying in the space and time. Finally, the effective action reads: 



S eff [A, A] «S +(a + c^) | A | 2 +^ | A | 4 . 
5.5 Time Dependence of the Linearized Equation 



(79) 



The linearized dependence requieres a careful examination of Gaussian fluctuations about 
the broken symmetry solution, and a simple low frequency expansion is obtained when 
A (T) << uj. In this case, there is necessary to expand Q (iqi) = r _1 (q = 0, iqj) — 
r _1 (q = 0, 0) in powers of uj after analytic continuation. From Eq. (70), we have: 



Q {Hi) = ^tanh 



1 



1 



2£k - iqi 2£ k J 



30) 



Under iqi — > a; + ie, and by using the relation (x ±ie) 1 = V (|) =F 2?r5 with P as the 
principal value, the frequency- dependent component becomes: 



tanh 



UJ 



2 72^(2^"^) 



+ i7r^tanh^^<S(2£ k -w). (81) 



The operator sum can be transformed into the integral through ^ k (• • •) = 

/ deN (ef) e~pl 2 E X l 2 (• • ■). Straightforward calculation leads to: Q (uj) = Q' (uj) +iQ" (uj), 

where: 

Aj( r \ roo cu£ 1/2 tanh ) 

Q' (uj) = ^M? / de- 

V 1 ; e F Jo 2{e-n)(2e-2n-uY 



^„ . , TriVfe^) , / uj \ Juj + 2u . „ , 



(82) 
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In the low frequency limit {u « T c ), an expansion of Q is possible for both the BCS and 
BEC regimes, where the condition uj «\ | is satisfied. Linearization of Q finally leads 
to: 

Q {uo) = dcu, (83) 

with 

^ tanh(/3A/2) np c N (fgOyg Q ,„ , 
d = JF2 + 1 7V2 ■ ( 84 ) 

k OE F 

The inverse Fourier transformation in the variables (x, t) allows to write the TDGL equa- 
tion: (q 2 — > —V 2 , uj — > —id/dt) 

V 2 d \ 

a-c— + b\A\ 2 -id—)A{x,t) = 0. (85) 

— III (_/ v J 

For weak coupling approximation, the cut off frequency is much larger than T c but con- 
siderably smaller than £p . The energy cut off is equivalent to truncating the Matsubara 
sums at a limiting number iVo = Qbcs/ZkT >> 1 [17 . The parameter a can be therefore 
rewritten as (for excitations near to the Fermi energy Ep): 

(3^1 Je^ {iuj n + e) {iuo n -e) (3 ^ n uo n y ' 



n>0 " v x v " ' v " ' ' n>0 

The order parameter vanishes at the critical temperature T = T c , therefore the following 
condition holds: 

where fines is, in general, the characteristic cut off frequency in the superfluid state. 
Then 

^11^1 1 T c 

The combination of Equations (71) and (88) leads to: 

a = N(E F )\n^. (89) 
For the non-linear coefficient b, we recall the expression (76) as: 



+ E 2 ) 2 ' 
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After integration over the upper-half complex plane in the variable s, we find: 



b K -N M 



— y 



Using (n + l/2)~ z = (2 Z — 1) ( (z), with £ (z) as the Riemmann-Zeta function, finally 
we obtain the quartic contribution term on TDGL: b = 7N (ep) ( (3) /8n 2 T 2 . In analo- 
gous way, for weak coupling interaction, the parameter c can be obtained from Eq. (70). 
The term which contains q is defined through the expansion of: J2k {^ u m — 60 1 {i^m + 
£k) 1 (1+5 1 / (^m + £k)) \ where (7 is defined in (73). For small values of the g-momentum, 
the relevant term in the last expression is those which contains the combination: 

2m" 1 (k • n) 2 (q 2 /2m) . 

Then c ~ T c Y2k (^m — £k) _1 (iw m + £k)~ 3 2m~ 1 k 2 cos 2 Evaluating the integrals in the 
k-space, the following result is obtained: c = 7( (3) N (ep) Ef/\2T 2 ti 2 . From Eq. (82), 
the main contributions in this limit are given in this by: Q'^ ~ uN(ef) /4£p, and ~ 
c^iV (ep) /8T C , (/j, rs £p), which define the coefficient d. 
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6. Calculation of the TDGL Coefficients for T <T r 



The previous analysis can be extended for the broken symmetry state and high-order 
expansions in the effective action S e ff [A, A]. This formulation permits to make physi- 
cally motivated approximations for all values of the coupling constants at low tempera- 
tures (T « T c ). We consider a uniform static saddle point A (x) = A which satisfies 
SS e ff/SA = 0. The relationship between the interaction parameter g and A can be 
obtained by using the Gaussian integration technique. The final result can be expressed 
as: l/g = ^ k tanh {(3E\ i /2) /2E^ where = a/££ + Aq (see Eq. 36). The regularization 
condition for the ultraviolet divergence in the gap equation is established by expressing 
the final result in terms of the s-wave scattering lenght a s , which describes the two-body 
attraction g in a low density system. The renormalized gap equation then reads (Cfr. Eq. 
37): 

m 1 tanh(/3£ k /2)- 



4ira. 



2Ev 



(90) 



To investigate the fluctuations around the saddle point A , we write A (x, r) = A + 
r) (x, r), and we recall the generic definition for the Green's propagator in Eq. (63) 



G- 1 = G 1 (Ao)+w+ + ^-, 



where Gq 1 (A ) is defined as: 



Go (A ) 



Q k (A ) T k (A ) 

^(A ) -s fc *(A ; 



(91) 



(92) 



where Q k (A ) = - (iu m + f k ) / (< + E%), and T k (A ) = A / (< + El). The effective 
action functional must be rewritten in terms of their mean field values A , A and their 
corresponding fluctuating fields 77, 77. Hence: 



S eff h r l) = Yl 



-1) 



Tr In Go 1 ( A ) + V ^^Tr (G ( A ) S) B 



n=l 



where £ is defined as: 



V q 

Vq 



(93) 



(94) 



The second-order combination in Tr[GSG'S'] can be expressed in terms of the spinors 
fields = [r]*,i]^ q \ and the correlation propagator IT through the combination rf q II q r) q , 
where IT is defined by 



Z~2k ^k^k+q ~ Ylik Gk^k+qJ 



(95) 
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It is useful to define the spinors compounds like: rj q = (X q + i9 q ) / y/2, with the properties 
A* = \- q and 9* = 9- q , which are usually identified as the phase and amplitude fluctu- 
ations, respectively. Under this transformation, the quadratic form is transformed like: 

n\%n, = (K K) fl) (£) , (96) 

and 

nf = \ (nf + nf + n- + nf ) , nf = 1 (nf - nf + nf - nf ) , 
nf = \ (nf ~ n f - nf + nf ) . 

The components n* J are defined in Eq.(95) in terms of the correlations functions — ^ k GkQk 
Fk^k+q and their complex conjugates. Our next task is to obtain approximate ex- 
pressions (second order in q and in the bosonic frequency oj), for each U.^ . 
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6.1 Fluctuations in the static Approximation: iv n = 0. 



We identify the following properties in the Green's tensor components at go 



The quadratic term in the static approximation for the fields A q and # q reads explicitly 
as: 

$Vfc (A* 0$ ( U * + U « nll ° ^ 2 ) fa) . (97) 



TlTl+ q (for A 



A n ). Hence: LI 1 , 1 



n 22 n 12 



V 



n 21 , and n q 2 



0: 
0. 







n 11 

q 



n 12 

q 



Or 



The coefficients in the diagonal form in Eq. (97) can be reduced as: 



12 



K ± n q 



^k +q gkCk+q ~ El T A 2 ) tanh (ggk/2) - E k fejk+q - ^ +q T A 2 ) tanh (/?£ k+q /2) 
^ 2£ k £ k+q (£ 2 - £ 2 +q ) 

For small values of the bosonic momentum q, (q 2 /m << min k {i? k }) the last expression 
becomes: LI 11 ± LI q 2 m c± (Ao, T, u) (q 2 /2m), where c± is formally defined in Appendix A. 
The particular case in which T = T c , (A = 0) the last results reproduce the expressions 
for a and c in Eqs. (71) and (75) respectively. In the opposite limit (i.e. T = 0), and 
for q = 0, the set of Equations (98) are reduced to: ITq 1 — LIq 2 = 0, which ensures that 
the phase mode for q is gapless, i.e., the Goldstone mode, and LTq 1 + LIq 2 = A 2 ] /2E^ L , 
which correspond to the results obtained by Engelbrecht et-al in Ref. [10]. 
The effective action Eq. (97) may be scripted as: 



s eff h V\ = S + (^q *) 



(99) 



where the matrix 1 



n q ) is diagonal in this approximation, and Sq 



Y2 q TrlnG 1 (A ). Let us define the generating functional Z for the bilinear form as: 



Z = Z Q / V[n, rf] exp 



ZoJJdet (tf q ). 



(100) 



where Zq = exp (—So), and So is defined in the textline after Eq. (99). Quantities of 
physical interest such as the free energy can now be directly computed by using T = 
— Tln£[l8j. Inserting the compounds of into the second term in the r.h.s of (100), 
we get: 



1 



T (A , fi, T) = To (A , fi, T) + - ^ In 



(9 



-l 



(101) 
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where JF = /3 _1 So, while the relations for 11^ are explicitly encoded in Eq. (98) and 
in Appendix A. In order to obtain the quadratic average fluctuations in the field 77, we 
consider the path integral representation in terms of the amplitude and phase components 
as follows: (notice that S does not depend on the elements K^ 1 ) 

\nZ = -S + J2^ Z <v Z^ = Z\Z^ (102) 
q 

where the functions Z x , Z 6 are respectively given by: 

Z\ = fv [A*, A q ] e-^)" 1 !^ 2 , Z% = Jv [0*,0 q ] e-^ 2 ^"l^l 2 . 
Therefore, there is straightforward to show that: 

or 

(I r) | 2 ) w=0 = 4 Y y - \ q (104) 

Y (^ 1 + n q 1 ) 2 -(n-) 2 

The numerical evaluation of this quantity as a function of the temperature will be dis- 
cussed elsewhere. An improved approximation for the order parameter is defined by 
adding the fluctuating term rj obtained in Eq. (104): A = A + a/(| f] | 2 )- This correction 
can be interpreted as a "pseudogap" the context of simple MFA theories. 
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6.2 Dynamics of Fluctuations: A Generic expression 

For the case in which the set of bosonic frequencies iv n is different of zero, the dynamical 
approximation is taken into consideration directly from the formula (99). The tensor K~ l 
might be separated into its "static" component, which has been collected into K~ l (See 
the textline right after Eq. 99), and its frequency-dependent part denoted as K~ l , where 
uj is the real frequency that appears after analytic continuation on iv, under iv — > u + iO + . 
Hence, K' 1 = K^ 1 + R-\ with 

rjr m-i ^ (l-A^)tanhQ3E k /2) 

[K -> 2E k (1-^/4^) ' ^ 

(K^y^-T tanh(/ ^ k/2) doe) 

(Kri2\-i ^ k tanh(/3E k /2) 

=t 2^ 4^(1-^/4^) ' (1 ° 7) 

with (JT, 2 , 1 ) -1 = — (If* 2 ) -1 . The expression (100) for the partition function Z might be 
modified as: 

\nZ = -S + J2 [ dujln l det ( K q)]> ( 108 ) 

where 

det (iC) = 



5- 1 + ni 1 + n-+(^ 1 )- 



5 



12 



+ n q -n q 



+ w 



1 

+ 4 



(A'i 2 )- 1 ! . (109) 



The last expression provides the core result in this report: The complete set of self- 
consistent equations that describes the Helmhontz free energy associated to Gaussian 
fluctuations in the field 77, via the relationship T (A , T, u) = In Z, for any tempera- 
ture between T and T c , and under the regularization condition in the BSC-BEC crossover 
scenario (Eq. 90). A similar expression to Eq. (104) is obtained for the full dynamics in 
the fluctuating field rj (u>, q): 



(hi 2 ) 



det {K-l) ' 



(110) 



where Tr{-} indicates the Trace over the embraced quantity. 
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Concluding Remarks 

In this paper we have shown how to formulate and derivate the quadratic Gaussian fluctu- 
ations associated to the perturbations in the order paramerter in a uniform Fermi gas, for 
any temperature below T c , by using path integral formulation. The average estimation 
for these fluctuations was obtained in exact form for the static case, while the partition 
function was drafted in a general way in terms of the BCS correlation propagators for the 
dynamical approach. 

The possibilities for a further work are listed in the following terms: 

> Numerical analysis for the final expression for (| 77 | 2 ), as a function of the tempera- 
ture and/or different values of the chemical potential //, and the collective excitations 
associated to the phase and amplitude modes in the BCS-BEC crossover region. 

> Perturbative analysis up to order four {O (| 77 | 4 )) in the effective action functional 

S eff- 

> TDGL treatment for trapped Fermi gas. 

> A long-term project might certainly be orientated in obtaining analogous expressions 
for a Fermi gas in the Feshbach resonance scenario |19j . 
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Appendix A: 

Definition of the "kinetic" coefficient c± in Eq. (98) 

The expansion in the second order in the wavector q leads to the expressions 2 : 



c± (A ,T,/i) = 



AEl 



with 



El 



+ 



(k ■ ny 



A 2 / A 2 
_1 + 10^2 ( l _ ±±° 

+ V El 



a , (k-n) 



mEh 



- 1 + 3 t>~ 



7'+ 



A 2 £ k (k • ii) 



AEi 



+ 



AmEv 



A 2 A 4 " 
1-4=^ + 6=^ 
El +b Et> 



AmE v 



and X = tanh (/3£/2), F = sech 2 (/3£/2). Notice that for A = 0, Z+ = Z", Z'+ 
and this results will reduce to those obtained by Randeira et-al. [8] . 



2 There is useful to consider the following results: (E = E^, £ = £k, E' = -Ek+q, £' = £k+ q ) 
E' 2 -E 2 ^Ax(l + x)E 2 , x = Q/4E 2 -Q 2 /16E 4 , Q = 2/i'£ + h' 2 , ti = £' - £, h! = (q 2 + 2q (k • n)) /2m. 
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Appendix B: 
Collective Modes 



By using the results in appendix A, there is possible to obtain the relationship for the 
collective excitations u (q) by solving det (K~^j = in Eq. (109). It is clear that for 
small values of momentum and frequency, 



-Ruj 2 + c. 



iBuo. 



with A = £ k A§X/2£ 3 , D = £ k (l - & 2 /E 2 )X/8E 3 , R = E k x / 8 ^j and B = 
J2 k £X/4E 3 , with reduced notation in the coefficients A,B,D,R, and c±,X defined in 
Appendix A. Therefore: 



-Ru 2 + c_ ( ±- 
2m 



bW = o. 



Neglecting terms of the form u 4 , q 4 and uj 2 q 2 , we get a first approach for the collective 
excitations modes: u> (q) = v s (T, A , /x) | q |, with the "sound velocity" given by: 



v s (T, A ,/x) 



Ac. 



2m (B 2 + AR) 



For arbitrary coupling, the integrals involved in A, B, etc, have to be evaluated numer- 
ically. In the weak-coupling limit, for all integrals peaked around the Fermi surface at 
T = 0. In this limit, the sound velocity takes the value: v s = v F /V3, (vf as the Fermi 
velocity) which corresponds to the well known result for the phase mode. 
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